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SUMHABT 



The ▼ariouB equatlone for eoDputlng profile drag "by 
the momentum method aire examined, and the errors arising 
from complete or partial neglect of compreBslhlllty ef- 
fects in the Jones equation (B . & M. Ho. 1688) and the 
Blcknell equatlen (ITACA Bep, Vlo. 667) are eTaluated. The 
Integrating method of Sllversteln and Katsoff (Jour.,AerOr 
Sciences, vol. 7, no. 7, May 1940) Is shown to he accurate 
over a vide range of Mach number and wake shapes. 



IITTBOSUCTIOH 

Iquatlons for the computation of the profile drag of 
a hody from an integration of the momentum defect In Its 
wake were developed independently hy A. Bets In 1925 (ref- 
erence 1) and hy B« N. Jones In 1936 (reference 2). Some 
of the assumptions made In the equations are not rigorously 
correct, hut It has been shown that the errors involved 
are small- (references 3 and 4). Both Jones and Bets ne- 
glected compressibility effects, And It was necessary to 
adapt their equations to compressible flow becnuse of the 
high velocities encountered In modern wind— tunnel and 
flight testing. This adaptation was made by Blcknell 
(reference 5), by Silversteln and Eatsoff .(reference 6), 
and by Wright In some unpublished notes. One of the pur-^ 
poses of this pauer Is to compare the results ' obtained 
frcm each of these methods and to ' Investigate the varia- 
tion caused by the differences In aesumptlbna and' pro— 
c e dur e • . 

■ ■ 

The use of any of the momentum drag equations cor- 
rected for compressible flow Involves a great many compu- 
tations because a rather complicated ■ expression must be 
calculated at a number of pelnts In the wake. Consider- 
able computation time may be ellmlnatdd by use of the 
integrating method outlined In reference 6. This method 
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evaluates the pr oport ionallt7 factor that exists between 
the drag and the average total. head Iobb In the wakn , a 
value which can he readily deliermlned from graphical or 
numerical Integration of point— hy— point measurements, hy 
use of an Integrating manometer (reference 7), or hy 
means of the averaging rake descrl'bed In reference 6. 
Certain assumptions, principally that of a predetermined 
wake shape, are necessary for the evaluation of the pro- 
portionality facter. It is the' purpose of this paper to 
determine experimentally the error entailed In the uen of 
this method over a wide range of Hach numbers and wake 
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The eq.uatlons of Bets and Jones differ because of the 
original assumptions, but there Is no significant differ- 
ence between the resulte (reference 4). Jones' method has 
been used in the further developments because of the 
greater simplicity in computation procedurp. 

The defect in mosentum caused by the drag of an air- 
foil is: 



assuming that the static preseures at stations C and 3 
(fig. 1) are equal, which is true if station 2 Is an in- 
finite distance downstream. (All symbols used In this 
report are defined in figure 1 and In the appendix.) As 
the measuremen.t of velocities at- great distances from a 
body is impractical, It Is necessary to convert equation 
(l) Into one with terms that are measursble^ In order to 
do this, Jones made the assumption that there Is -no mitr- 
ing in the tubes of flow between plane 1, a measurement 
plane close behind the airfoil, and plane 3.. This assump- 
tion permits the application of Bernoulli's equation be- 
tween the two planes. 



shapes . 




(1) 



Then 



T^i-Ps ' P 




= Hj-p^ 



(2) 



and equation (1) reduces to 




(3) 
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or 



-.3 r -^o-Pi-dloT^Hi) ' 



As the' Jones eqliat Ion .,1a appllca'ble only to Ineom— 
pjresB 11)10 flov, fh^ nnxt r ^If^inemnnt was to adapt It to 
compresBl'ble flow. Blclczlt^ll (r^f^rencn 5) ehovs how 
Bernoulli's equation f or 'compresa Ihle flow is applied to 
determine' ^ and how the ohange In density caused hy 'the 
atfttle pressure diifferenca Iretween stations 0 and 1 may ' 
"be taken 'into acoouat hy thd' asfsumption of adlabatlo 
variations. With these modifications, Jones' equation 
he comes : 




Blcknell did not'conalder the increaaein temperature 
In the wake land the cohaequent density changea-, an effect 
which is appreclahle at high Mach numhers. This tempera- 
ture rise was taken' into account by &ay H, Wright in some 
unpublished work •and hy Silversteln and Katzoff in refer- 
ence 6. The same aa'atrmpt lone are made In hoth develop- 
ments, and the 'methods differ only in procedure. 

On the asaumptlon that the flow in eac^h streamline 
tube In the- Wake betwee-n places '1- and 2 'Is isentroplc'; tho 
tomperature rise in each lamina ;1b oqueil to th»' ^neiigy- dif- 
ference betweon the work exifendeld ■ pnr unit time Iz^. over- 
coming the drag and the klnet Ic ..eno.rgy produced in' the wake. 
' The aaaumption of 'Ise'ntropic flow la an&logous to the as— 
- a'umptlon fflade_ in the development of the inQompreaslble flow 
eqiiation that there' is no mixing in 'the tuhe^ 'of fl'ow be- 
tween .'planers 1 and 2- cLnd , consequently,, that, the t otal head 
rnmaln'e' constant (r of e']!-eno<! 9). ' ^per.imental verification 
of the 'isentrop'lc nature of the wake flow iis. given. in ref- 
erentJW 1-0-, "the resii'llii/'iof which iBh,ow .^th^t.; ft c.on8tant stag- 
nation' temperature exl.stB aor OB ^ '.th^e. yake , .tls.ihg t.he rela- 
tions' given in reference 6, the 'f o;Ll.oK;Lng equation, can'be 
developed' for the profile dr.ag .jcoe'f^f icient . (See ap^^hdix.) 



^ (Ho-Po)(i+ni) 



l+0.2.0SMa 



dy 



(6) 



■Vrlglit developed a almllar equation , tut Introdneed the 
.faotor 0 to define the relation '^■pa/po " ^"^(l-^a/^lo^ 

where 



0 « 



"v 1+Q.202M2 • 

1 - qa/qp 
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-Pi 



(7) 



0 Is. primarily a funotion of Maoh num'ber , and Beoondarlly 
a funotion of (l— qa/qo)* variation i of 0 with these 
two factors is shown in figure 2r 

I 

. Ta'ble I and figure 3 compare the profile drag.looef— 
f loients as computed from point-'l)y— point methods -from wake 
measurements on a 66,2-420 airfoil (reference 7} through- 
out .a^n extenelTe' Hach num'ber range* As expected, the 
iTones equation for incompressible flow gives the highest 
drc^g coefflolient at any Maoh number. Blcknell's equation 
Is Intermedl&te between compressible and .Incompress Ible ■ 
flow becausis'the temperature rl8^.|ln the twake was neglected. 
The two corrected methods of. computing .the drajs in oom-^ 
ptesslble fl.ow (equations (6) and . (7)T 'filTO results which 
differ only because of inaoouraoieff in -'plotting and com- 
puting. ■ '. , . 

A oouparisei^ of the drag ooefflcients at a high Mach 
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number glTRB an indication of thn magnitudea of thr> various 
correetions. .At- ..H ■= 0..6^ the, percentage differi^nce "be- 
tween Jones' and Bicknell's equations is 7 percnnt . Indi- 
cating the effect of eonpr i>s s 1*6 il ity , the inclusion of the 
l+r\ factor, and th^ density changos due t-e static pressure 
d if f erencPB . The effect of tho t emperaturfl rlsf in the 
wake on the density causes a further correction of 6 per- 
cent, as shown .the percentage differen-ce bistween 
Bicknell's and the Si Ivers te in— Kat spf f or Vri^ht equation. 
The total value of the corrections to the Jones equation 
at this Mach number is 1? percent. Assuming, that equations 
(6) aivd (?) give fully corrected reetults , the error from 
the Jones equation will exceed 2 percent if it is used abov«* 
a Mach number of 0.1,' and the grror from Bicknell's equation 
will exceed 2 percent If used above a Hach number of 0.25. 
The error In the Joa«>B equation increases roughly as M^-^; 
that of the Bicknell equation incrpases appr oxiuately In 
proportion to M. 

It should be noted that the drag results determined In 
NACA high-speed wind tunnels by the momentum method have 
included all of thf> compressibility and temperature correc- 
tions outlined above. 



IKTSGBATIua METHODS 



Equations (6) and (7) are somewhat unsa't isf act ory for 
vind-tunhel or flight test work because of fhe lengthy com- 
putations involved; the equations must be evaluated for a 
number of points throughout the wake, and the results 
graphically or numerically integrated. As developed in 
reference 6, a simplification results from the feet that 
the profile drag coefficient is roughly proportional to 
the average total— head loss In the wake. 



or 



(H — H ) ■ 
**d« "Pi - — (incompresslblei flow) (8b) 

O H — T) 



Cd^ * ^c r — fi — (c ompr ess ibile" .f low) 



^0 ° a H„--D„ 
0—0 



(8c) 



The value of (Hq-HO may be determined from a graph- 
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leal or numerical arerage of point— "by—point meaaurements , 
d-lreotly "by use of an integrating manometer (reference 8), 
or "by the averaging rake deacri'bed in reference 6. 1b 
a proportionality factor for IneompresBible flow and is a' 

Px-Po f^o"^i5max 
function of . • . and vake shjape. It is 

Sq-Po Hq-Po 
evaluated in reference 6 on the aasumpt ion that the wake con- 
forms to a cosine— squared curre. Xquations (3) and (61)) 
are evaluated In terms of this assumption and equated; then 
Tj^ Is determined for a range of values of Pi—Po/^o'Po 
and (Ho-H^)Bax/^o~Po* variation of with these 

two factors in the rangs normally encountered is shown in 
figure 4. 

To correct the Integrating method for compress 11)11 ity 
and temperature effects, Sllversteln and Eatsoff (reference 
6) determined Tq, the proportionality factor in compree— 
slhle flow, hy evaluating equations (6) and (Sc) In terms 
of the assumption of a cos ine— squared wake form, equating 
and solving for Vq. Values of I'p/l'i so determined are 
presented in tahle I or reference 6 as a function of 
"(Ho-Pi)max Pi-Po 



=o-Pn H 



and Mach number 



o 



In using this method In wind— tunnel testing, It has 
heen found convenient to determine 7^ first and then 
modify It by the ratio ^o/^i ^ °^ Mach numher of the 

test velocity. Thus, 

'o 



' 'o ? ' T ^^ ^ ^ (8c) 



'l !- ^ ° =o-Po 



(9) 



The variation of 'q/'i vlth the three variables involved 
is given in figure 6. 

It is evident from equation (9) that any error in- 
volved In the assumption of a c os Ine— squared wake form is 
primarily a function of and secondarily a function 

of Fq/Tj^. In order to eliminate any error in T^, the 

following equation may he used since 'o/^i shows the 



7 



relation betveen oofflpress 11)16 and InoonpreBslble flovt 

I Drag as deterialne4"b3r the " "j " " (io) 

^0 " j_Jonea equation (ec^uatios (3))J 

the aoouraoj of the Integrating method was inyeeti'- 
gated from the wake meaflurementB o'btalned In reference 7« 
The profile drag ooeff iolente were determined hy point— 
hy—point oomputat ions and the integrating method from the 
wake BtirTeys vhioh were made through ^ Maoh number range- 
of O.iSB to 0.668. Taible IX compared the results and. . 
shows that the ooeffioients determined by equation (9) are 
the same as those obtained in equation (10). In other 
words, eontalns very little .errQr for. the wake form 

enoountered in the teste, Compariesn of the results of 
equatione (10) and (6) show that the error in 'q/'i Is 
also negligible. At all speeds Tieiew the critioal, the 
greatest difference between the integrating method and the 
point— by—point computation of equation (6) was 1.7 percent* 

Wakes may vary considerably in shape} so several dif- 
ferent forms were assumed, to determined the effect of ex- 
treme departures from the ooe ine—Bquared assumption. The 
profile drag coefficient for an -unBymmetr leal , a triangular, 
a rectangular, and a true oob Ine— squared wake form (fig. 6) 
were computed by the point— by— point methods and by the In- 
tegrating method ueing 7j^^ and Tq/J^, The uneymmetrloal 

wake le the type found behind a wing— nacelle Junction and 1b 
due te the combination of the wakee of each. Table III 
Shows that the actual shape has little effect upon the 
error in the profile drag coefficient as cemputed from the 
COB ine— squared asBumptien. The percentage errsre between 
the integrating and the point-i-by— point computation methods 
show the error In , 'q/'ii and the produot ef ' the two. 

There appears to be no consistency as the wake form departs 
more and more from the cos Ine— squared shape; the major por- 
tion of the error is prebably due to Inaqottraoies in plat- 
ting and csmputlng. 



GOHCLUSIOHS 



1. The Jones equation has less than a 2— percent err«r 
due te compressibility effects, up to. a Maeh number ef C.l. 
Thi;8 error increases rapidly With Mach nixinber , jfoughly in 
proportion te k*-". " > 
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2. The Blcknell equat Ion. hae less than a 3— peTcen-t. 
error up to a Mach number of 0'.25 due to compress Ih lilt/ 
effects. This err or ,■ Increases approximately In proportion 
to- Mach numher . , 

f 

I 

3. The integrating method of Sllvereteln and Eatsoff 
gives results which agree within 2 percent with the.re~ 
suits of the integration of point— b/— point computations 

. (corrected -f or compreBs iblllty effects) at least up to the 
CFi:t leal ' Mach number for normal wing wakea. The . addlt lonal 
error 'due to extreme variation in wake shape is less than 
1 . 5' percent . 

Ames Aer onaut'ical' Laboratory, . 

national Advisory Committee for Aeronautics, 
■ ^ Moffett rie.l4," Calif . .'■ 



APPENDIX 

.'t.t 't':- ' .«•••••..''■■ 

Definitions of symbols used in this r.eportrf , 

-.■I 

q. . dynamic pressure — l/2p7' ■ H— p 
p static pressure 

* . -4.1 ■ ■ ■ 

■ 

H ^ tqtal prcissure 

V velocity 

.1 ■ "• 

p dens Ity 

* ■ ' ■ 
T absolute temperature, 

dS elemental area- 

. ■ ■ ■ .-r .- ■ 

y vertical distance from wake center 
w width of wake . 
c airfoil chord 

Y* ' ratio of specific :heat at constant pressure to\ihe 
' apqelf,lc heat at constant volume 
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H Mach number, the -ratio of . the etream yeloclty to the 
local Telocity of sound 

7 proportionality factorv The euhacript i -denotes - 
Incompresslhle flow and c, compr esslhle flow -- 

e^^ section profile drag coefficient . 

1+11 a coopressihlllty- correct Ion factor 

'. ' q. » = 1 eompresslhle flow 

where 1+ti = 1 + ^ + + 

.'The subscripts o« i', and a denote the location 
of the measured or computed values t o denotes the free 
stream; i denotes the measurement plane, close behind 
the airfoil; a* denotes the value at the plane located 
an Infinite distance downstream from the airfoil. 



The derivation, of .equat ion (6), according to the 
relations set up In reference 6,' is as follows: 

■ 4 - P£Vs(Vo-Va) dya (a^ 

|r £5l5(i ,Is)dya . (b) 

po^d ^ . ^0 ^ 



"d. 



w 

Tor continuity: PiV^dy^ => pa'^a^^s 



°-o ... o / PoVq \^ 



The difference, betveen the work expanded per unit 
time an^d the kinetic energy produced In t'he wake Is: 

PaVaVo(Vo-V8)dya - a q-^ a)^ 9 - |paTa (▼o'-'^a'')dya 

If this difference remains in the lamina as heat, the tern 
perai.ture rise 'Is;'- . ' . ' . 
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where 



T abaolute temperature; fteg. C 

w - i. ■ 

e specif lo heat at eonetant preB&u-re 

Suhstltutlon of the numerical value of o. and dlylelon 
Toy a 0.000231 V^.^* leads to ' • ^ 

Ta-T s / Va' \ 
. .2 = -CSOS (.1 1 - (e) 

Since p; = Pa. » -'-1 + 0.302 (l - I^," ) " ' (f ) 

Pa ti ^ a ^ *o' 



^o y^Pa (H -p )(1+Tia) y Pa (Hq-Pq) (l+rja) 

0 0 

Solving these tvo equations simultaneously for po/pa 

p 1 + 0.202 m' 

: (h) 

pa a/^i-PoN 



1+ 0.202 H-{ ; ( ^ ) 

VHg-p^y \l+tia y 

wher.e M a 7q/Tq^ and is the.Maeh number. 

li ,-/^ (gi-Pi) (1+Tlo )' pjTi ^ /pi (Hi^Pi)(H-tio) 

■Po Pa P.O. . 
Since the flow between planes 1 and 2 is isentroplc: 
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/Pi >Y /PiN.Y 

at; A 



Pi _ 
Pa ^Pa ^■ 



(k) 



Substituting equations (k) and. (f) In. eq.ua^t.ion (j) and 
(J) into (i):. 



0^0 ' -A^^o^ (Ho-Po)(l+^i) V 1 + 0.202Ma 

■ ■ 1 I 



Pn)(.l+ni) 



''(Hi-p„)(l+T,„) 



1 + 0.302 MS y (Ho-PoXi+tia) 

Wright used the relation ■/ps/po ■ 1 - (l-QLa/q.©) 



dy (6) 



where 



0 = 



1 -J>Jp^ 



k ^ 0 .202 
y 1+0. 



202M» qa/q, 



202M" 



1 - <la/tlo 
to obtain equation (?). 



1 - la/4, 



2 r /Pi\:^ (Hj^ 
° Jw y^PpV .. (Ho'Po 



-Pi ) (l+mp) 
).(l+ri,). 



_afi-ii:££^V 

\ . Hq-Po l+TJa-^ 



1+Tl, 



dy 
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TABLE 


I 




V8 M as eofflputed "by the varloue momentum drag 
0 methods from tests on a 66,2—420 airfoil 


number 


J one B ' 

eq.uat 1 on 
\o ) 


V J nVn all la 

equat Ion 


Q 4 1 w A w a ^ A 4 n 

3 1 XV Sir 8 W 6 1 

Eat feoff 
equabiou ^D^ 


HT Igllli B 

eauat ion 

\r ) 


0. 188 


0 , 00568 


0 . 00564 


0. 00556 


0 . 00556 


.271 


.00637 


.C0628 


,00616 


.00618 


.328 


.00629 


.00622 


.00612 


.00613 


.3 83 


.00637 


.00623 


.00607 


.00603 


.515 


.00698 


. 00659 


,00633 


.00637 


.594 


.00756 


. 00703 


.00662 


.00656 


.63 8 


.00888 


.00838 


,00773 


.00771 


. .658 

1 


.01218 


. 01152 


. 01066 


,01031 
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TASLS II 






cj tb M to determine 

0-0 


the error 


In ■ 


?! and Pr/Fi 


Mach 
number 


Integrating 
method x 
X Pg/Fi 


equat 1 on 

V /in 
X Fg/Fi 




r** ■ 

1 

1 

SllTPrstein— 
1 Eat z off 
1 equation (6) 




equation (9) 


equat 1 on 


(10) 1 ■ 


0.188 


0.0055 9 


0. 00560 




0.00556 


.271 


.00621 


. 0C618 




.00C16 


.R38 


.00604 


. 00602 




.00612 


.?83 


.00599 


.00600 




. 00607 


.515 


.00632 


.00629 




.00633 


.594 


.00662 


. 00o5 9 




.00662 


.63a 


. 0075S 


. 00759 




.00773 


.658 


.01010 


.01034 




.01066 
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TABLZ III 








0 


wake form to 


determine the error In and 'e^'l 




Vake 
form 


SllverBtein— 
Eats off 
equation (6) 


Jones ' 
equat Ion 

X 'cAl 
equat Ion 
(10) 


Percent 
error 


Integrating 

method 

X Tc/Ti 

equation (9) 


Percent 
error 




cos^ 6 


0. 00698 


0.00899 


+0.1 


0.00904 


+ 0.7 




Trian- 
gular 


.00912 


.00902 


-1.1 


.00904 
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Wing- 
nacelle 


.01314 


.01295 


-1.4 


.1308 






Etectan— 
gular 


.01786 


' .01773 


-.7 


,01807 


+ 1.2 



InconiproBsl1)le flow 



Piano 0 



CoDiprosBl'blo flotr 



'o'^o'^a'Po'^^o 
gp-Po 



%'Pl'^l 
^o'Po 



Hi.Pq.T2 
''o.Po 

^aTHi-Po 




Piano 1 



Piano 2 



7 = Vortical dlBtance 
from Tniko contor 



Tl.Pl 
^1 1+Tli 



ifS'Pg 



Figaro 1.- Boflnltlon of toxms uaod In tho momjntam oqoatlons 



IPig. 3 



O JonoB ' ognaatlon 
^Blcknoll'B " 
a Wright 'B n 
+ SatBoff and SllvorBtoln'B 
oqaatlon 
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.3 .4 

Uach. numbor 



Vlgure 3.- ys. Uach namber by the varlciaB polnt-'byopulnt methodB. 



(Mmasut-m with '4<o') 




ng. 4 



(Ho-Hi)muc 

Figure 4.- ?! (Ho-Bi.)aax/Bo--Po for oonataat Talusa of Pi-Pe/Bo-Po* 




.760 




Pig. 6 




VeOeo Trldth - Inchoa 
figure 6.- Wako forma. 



